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Abstract 

The ARGO-YBJ experiment is a full-coverage air shower detector located at the Yangbajing 
Cosmic Ray Observatory (Tibet, People’s Republic of China, 4300 m a.s.l.). The high altitude, 
combined with the full-coverage technique, allows the detection of extensive air showers in a wide 
energy range and offer the possibility of measuring the cosmic ray proton plus helium spectrum 
down to the TeV region, where direct balloon/space-borne measurements are available. The de¬ 
tector has been in stable data taking in its full configuration from November 2007 to February 
2013. In this paper the measurement of the cosmic ray proton plus helium energy spectrum is 
presented in the region 3 — 300 TeV by analyzing the full collected data sample. The resulting 
spectral index is 7 = —2.64 ±0.01. These results demonstrate the possibility of performing an 
accurate measurement of the spectrum of light elements with a ground based air shower detector. 


* Corresponding author: stefanomaria.mari@uniroma3.it 
1 Corresponding author: paolo.montini@roma3.infn.it 
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I. INTRODUCTION 


Cosmic rays are ionized nuclei reaching the Earth from outside the solar system. Many 
experimental efforts have been devoted to the study of cosmic ray properties. In the last 
decades many experiments were focused on the identihcation of cosmic ray sources and on 
the understanding of their acceleration and propagation mechanisms. Despite a very large 
amount of data collected so far, the origin and propagation of cosmic rays are still under 
discussion. Supernova remnants (SNRs) are commonly identihed as the source of galactic 
cosmic rays since they could provide the amount of energy needed in order to accelerate 
particles up to the highest energies in the Galaxy. The measurement of the diffuse gamma- 
ray radiation in the energy range 1 — 100 GeV supports these hypotheses on the origin and 
propagation of cosmic rays [I]. Moreover the TeV gamma-ray emission from SNRs, detected 
by ground-based experiment, can be related to the acceleration of particles up to ~ 100 TeV 
A very detailed measurement of the energy spectrum and composition of primary 
cosmic rays will lead to a deeper knowledge of the acceleration and propagation mechanisms. 
Since the energy spectrum spans a huge energy interval, experiments dedicated to the study 
of cosmic ray properties are essentially divided into two broad classes. Direct experiments 
operating on satellites or balloons are able to measure the energy spectrum and the isotopic 
composition of cosmic rays on top of the atmosphere. Due to their reduced detector active 
surface and the limited exposure time the maximum detectable energy is limited up to few 
TeV. New generation instruments, capable of long balloon flights, have extended the energy 
measurements up to ~ 100 TeV. All the information concerning cosmic rays above 100 TeV 
is provided by ground-based air shower experiments. Air shower experiments are able to 
observe the cascade of particles produced by the interaction between cosmic rays and the 
Earth’s atmosphere. Ground based experiments detect extensive air showers produced by 
primaries with energies up to 10^° eV, however they do not allow an easy determination of 
the abundances of individual elements and the measurement of the composition is therefore 
limited only to the main elemental groups. Moreover, due to a lack of a model-independent 
energy calibration, the determination of the primary energy relies on the hadronic interaction 
model used in the description of the shower’s development. 

The ARGO-YBJ experiment is a high-altitude full-coverage air shower detector which was 
in full and stable data taking from November 2007 up to February 2013. As described in 
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section [TT| the detector is equipped with a digital and an analog readout systems working 
independently in order to study the cosmic ray properties in the energy range 1 — 10^ TeV, 
which is one of the main physics goals of the ARGO-YBJ experiment. The high space- 
time resolution of the digital readout system allows the detection of showers produced by 
primaries down to few TeV, where balloon-borne measurements are available. The analog 
readout system was designed and built in order to detect showers in a very wide range 
of particle density at ground level and to explore the cosmic ray spectrum up to the PeV 
region. In 2012 a hrst measurement of the cosmic ray proton plus helium (light component) 
spectrum obtained by analyzing a small sample collected during the hrst period of data 
taking with the detector in its full conhguration (by using the digital readout information 
only) has been presented [1]. 

In this paper we report the analysis of the full data sample collected by the ARGO-YBJ 
experiment in the period from January 2008 to December 2012 and the measurement of 
the light component energy spectrum of cosmic rays in the energy range 3 — 300 TeV by 
applying an unfolding procedure based on the bayesian probabilities. The analysis of the 
analog readout data and the corresponding cosmic ray spectrum up to the PeV energy region 
is in progress and will be addressed in a future paper. 


II. THE ARGO-YBJ EXPERIMENT 

The ARGO-YBJ experiment (Yangbajing Gosmic Ray Observatory, Tibet, P.R. Ghina. 
4300 m a.s.l.) is a full-coverage detector made of a single layer of Resistive Plate Ghambers 
(RPGs) with ~ 93% active area 13 E], surrounded by a partially instrumented guard ring 
designed to improve the event reconstruction. The detector is made of 1836 RPGs, arranged 
in 153 clusters each made of 12 chambers. The digital readout consists of 18360 pads each 
segmented in 8 strips. A dedicated procedure was implemented to calibrate the detector in 
order to achieve high pointing accuracy [7] . The angular and core reconstruction resolution 
are respectively 0.4° and 5 m for events with at least 500 hred pads m- The installation of 
the central carpet was completed in June 2006. The guard ring was completed during spring 
2007 and connected to the data acquisition system [10] in November 2007. A simple trigger 
logic based on the coincidence between the pad signals was implemented. The detector 
has been in stable data taking in its full conhguration for more than hve years with a 
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trigger threshold Npad = 20, corresponding to a trigger rate of about 3.6 kHz and a dead 
time of 4%. The high granularity and time resolution of the detector provide a detailed 
three-dimensional reconstruction of the shower front. The high altitude location and the 
segmentation of the experiment offer the possibility to detect showers produced by charged 
cosmic rays with energies down to few TeV. The digital readout of the pad system allows 
reconstruction of showers with a particle density at ground level up to about 23 particles/m^, 
which correspond to primaries up to a few hundreds of TeV. In order to extend the detector 
operating range and investigate energies up to the PeV region each RPC has been equipped 
with two large size electrodes called Big Pads m- Each Big Pad provide a signal whose 
amplitude is proportional to the number of particles impinging the detector surface. The 
analog readout system allows a detailed measurement of showers with particle density at 
ground up to more than 10^ particles/m^. 


III. DATA ANALYSIS 


A. Unfolding of the cosmic ray energy spectrnm 


As widely described in [12], the determination of the cosmic ray spectrum starting 
from the measured space-time distribution of charged particles at ground level is a classical 
unfolding problem that can be dealt by using the bayesian technique [T3|. In this framework 
the detector response is represented by the probability P{Mj\Ei) of measuring a multiplicity 
Mj due to a shower produced by a primary of energy Ej. The estimated number of events in 
a certain energy bin Ei is therefore related to the number of events measured in a multiplicity 
bin Mj by the equation 

N{E,) ^J2N{M^)P{Ei\Mj) (1) 

j 

where rjij is constructed by using the Bayes theorem 


P{E,\M,) 


P{M,\E,)P{E,) 

EkP{M,\E,)P{Eky 


( 2 ) 


The values of the probability P{Mj\Ei) are evaluated by means of a Monte Carlo simulation 
of the development of the shower and of the detector response. The quantity P{Ei\Mj) 
represents the probability that a shower detected with multiplicity Mj has been produced 
by a primary of energy E^. The values of P{Mj\Ei) are evaluated by means of an iterative 
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procedure starting from a prior value in which in the n-th step P^'^\Ei) is replaced 

by the updated value 


P^^^^\Ei) 


N^^\Ei) 


(3) 


where N^'^\Ei) is evaluated in the u-th step according to eq. ^ As initial prior P^^^Eij ~ 
E~‘^'^ was chosen, the effect of using different prior distributions has been evaluated as 
negligible. The iterative procedure ends when the variation of all N{Ei) in two consecutive 
steps are evaluated as negligible, namely less than 0.1 %. Typically the convergence is 
reached after 3 iterations. 


B. Air shower and detector simulations 

The development of the shower in the Earth’s atmosphere has been simulated by using 
the CORSIKA (v. 6980) code [H]. The electromagnetic component are described by the 
EGS4 code [l5l[T6], while the high energy hadronic interactions are reproduced by QGSJET- 
11.03 model [171118]. Low energy hadronic interactions are described by the FLUKA package 
[T9l|20]. Showers produced by Protons, Helium, GNO nuclei and Iron have been generated 
with a spectral index 7 = —1 in the energy range (0.316 — 3.16 x 10'^) TeV. About 5 x 10^ 
showers have been generated in the zenith angle range 0-45 degrees and azimuth angle range 
0-360 degrees. Showers were sampled at the Yangbajing altitude and the shower core was 
randomly distributed over an area of (250 x 250) m^ centered on the detector. The resulting 
GORSIKA showers have been processed by a GEANT3 121 ] based code in order to reproduce 
the detector response, including the effects of time resolution, RPG efficiency, trigger logic, 
accidental background produced by each pad and electronic noise. 


C. Event selection 

The ARGO-YBJ experiment was in stable data taking in its full conhguration for more 
than hve years: more than 5 x 10^^ events have been recorded and reconstructed. Several 
tools have been implemented in order to monitor the detector operation and reconstruction 
quality. The detector control system (DGS) [22] continuously monitors the RPG current, the 
high voltage distribution, the gas mixture and the environmental conditions (temperature. 
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pressure, humidity). In this work the analysis of events collected during the period 2008- 
2012 is presented. Data and simulated events have been selected according to a multi-step 
procedure in order to obtain high quality events and to ensure a reliable and unbiased 
evaluation of the bayesian probabilities. The hrst step concerns the run selection: in order 
to obtain a sample of high-quality runs, the working condition of the detector and the 
quality of the reconstruction procedure have been analyzed by using the criteria described 
below. 

• At least 128 clusters out of 130 must be active and connected to the DAQ and trigger 
systems. This criterium selects runs taken with almost the whole apparatus in data 
taking, discarding the runs that can bias the analysis because of the switched-off 
clusters. 

• Only runs with a duration T ^ 1800 s have been considered. The runs with a short 
duration are generally produced when a problem in the apparatus occurs. These runs 
have been removed from the analysis. 

• The value of the trigger rate for each run must stay within the range 3.2 — 3.7 kHz. 
A trigger rate outside this range usually indicates that the detector was not operated 
standardly. These runs have been discarded. 

• To monitor the quality of the event reconstruction the mean value of the unnormalised 

obtained by htting the shower front must be less than 135 ns^ (see £gure[^. Nearly 
all runs that have > 135 ns^ encountered some sort of problems. 

In hgure [T] the distribution of the trigger rate and the of the reconstruction procedure 
are reported. The procedure described above selects a data sample of about 3 x 10^^ events, 
corresponding to a live time of about 24000 hours. 

The following selection criteria (hducial cuts) have been applied to both Monte Carlo and 
experimental data in order to improve the quality of the reconstruction and to obtain the 
best estimation of the bayesian probabilities. 

• Only events with reconstructed zenith angles 'dn ^ 35° have been considered. The 
resulting solid angle D is about 1.13 sr. 


7 




FIG. 1. Distribution of the trigger rate (top) and of the unnormalised (bottom) of all the runs 
collected by the ARGO-YBJ experiment (black lines). The resulting 2008-2012 sample selected 


according to the criteria described in section III is also reported (dashed red lines) 


• The measured shower multiplicity M had to be in the range 150 ^ M ^ 5 x 10^. This 
selection cut was introduced in order to reduce bias effects in the estimation of the 
bayesian probabilities that are mainly located at the edges of the simulated energy 
range. Moreover the highest multiplicity cut avoid saturation effects of the digital 
readout system. 

• The cluster with the highest multiplicity had to be contained within an area of about 
40 X 40 m^ centered on the detector. This cut was applied in order to reject events 
with their true shower core position located outside the detector surface. 

In order to select showers induced by proton and helium nuclei the following criterium has 
been used. 

• Density cut: the average particle density (pm) measured by the central area (20 inner 
clusters) of the detector must be higher than the particle density {pout) measured by 
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FIG. 2. Distribution of reconstructed core positions of showers selected by applying the criteria 


described in section III C 


The boxes represent the clusters layout. 



FIG. 3. Energy distribution of all Monte Carlo events (black) and of those surviving the hducial 


cuts (blue) and the density cut (green and red) described in section III C according to the Horandel 
model 


the outermost area (42 outer clusters): {pin > 1.2bpout)- This selection criteria based 
on the lateral particle distribution was introduced in order to discard events produced 
by nuclei heavier than helium. In fact, in showers induced by heavy primaries the 
lateral distribution is wider than in light-induced ones. By applying this criterion 
on events with the core located in a narrow area around the detector center, showers 
mainly produced by light primaries have been selected. The contamination of elements 


heavier than helium does not exceed few %, as discussed in section IV A 4 
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FIG. 4. The light component spectrum measured by the ARGO-YBJ experiment by using data 
taken in each year of the period 2008-2012 and the full 2008-2012 data sample. The power-law 
fit of each spectrum is also reported (red lines). 

In figure the coordinates of the reconstructed core position of the events surviving the 
selection criteria described above are reported. The plot shows that the contribution of 
events located outside an area of 40 x 40 is negligible. In hgurej^the event rate obtained 
by using the Horandel model for input spectra and isotopic composition [2^ and surviving 
the selection criteria described above is reported as a function of energy for both proton 
plus helium (light component) and heavier elements (heavy component). The plot shows 
that the selected sample is essentially made of light nuclei. 


IV. THE LIGHT COMPONENT SPECTRUM 


The analysis was performed on the sample selected by the criteria described in section 


TTTj Simulated events have been sorted in 16 multiplicity bins and 13 energy bins in order 
to minimize the statistical error and to reduce bin migration effects. The Monte Carlo 
data sample was analyzed in order to evaluate the probability distribution P{M\E) and the 
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energy resolution which turns out to be about 10% for energies below 10 TeV and of the 
order of 5% at energies of about 100 TeV. The multiplicity distribution extracted from data 


has been unfolded according to the procedure described in section [III A[ Results are reported 
in hgurej^for each year of data taking and also for the full sample. In order to investigate 
the stability of the detector over a long period the analysis was performed separately on the 
data samples collected during each solar year in the period 2008-2012. The values of the 
proton plus helium flux measured at 50 TeV are reported in table |T| A power-law £t has been 
performed on the measured spectrum of each year and of the full data sample, the resulting 
spectral indices are reported in figure Both the spectral indices and the flux values are 
in very good agreement between them, demonstrating the long-period reliability and the 
stability of the detector. The spectral index 7 = —2.64T0.01, obtained by analyzing the full 
data sample, is in good agreement with the one measured by using a smaller data sample 
collected in the hrst months of 2008 | 1 ] which was not corrected by the contamination from 


heavier nuclei (see section IVA4). 

In table and hgure the flux obtained by analyzing the full data sample is reported. 
The spectrum covers a wide energy range, spanning about two orders of magnitude and is 
in excellent agreement with the previous ARGO-YBJ measurement. Statistical errors are 
of the order of l%o, more than 10 ^ events have been selected in the highest energy region, 
while at the lowest energies more than 10^ events have been selected. Systematic errors 
are discussed in the next section. The ARGO-YBJ data are in good agreement with the 
GREAM proton plus helium spectrum [21]. At energies around 10 TeV and 50 TeV the 
fluxes differ by about 10% and 20% respectively. This means that the absolute energy scale 
difference of the two experiments is within 4% and 6 %. The uncertainty on the absolute 
energy scale has been evaluated by exploiting the Moon shadow tool at a level of 10% for 
energies below 30 TeV [ 8 ]. At present the ARGO-YBJ experiment is the only ground-based 
detector able to investigate the cosmic ray energy spectrum in this energy region. 


A. Systematic uncertainties 

A study of possible systematic effects has been performed. Four main sources of system¬ 
atic uncertainties on the flux measurement have been considered in this work: variation of 
the selection cuts, reliability of the detector simulation, different interaction models, con- 
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FIG. 5. Spectral indices of the power-law fit of the light component spectrum measured by an¬ 
alyzing the data sample collected in the period 2008-2012. The spectral index obtained in a 
previous analysis of the ARGO-YBJ data is shown as 2008* [3j . The error bars represent the 
total uncertainty. 

TABLE I. Proton plus helium flux measured at 5.0 x 10'^ GeV. 



Year 

Flux ± tot. error [m ^s ^sr 

2008 

(4.53 ±0.28) X 10-9 

2009 

(4.54 ±0.28) X 10-9 

2010 

(4.54 ±0.28) X 10-9 

2011 

(4.50 ±0.27) X 10-9 

2012 

(4.36 ±0.27) X 10-9 


tamination by heavy elements. 


1. Selection criteria 

The fiducial selection criteria have been fine tuned in order to obtain an unbiased evalu¬ 
ation of the bayesian probabilities, leading to the best estimation of the cosmic ray proton 
plus helium energy spectrum. A possible source of systematic error is related to the values 
of the hducial cuts on observables used in the event selection procedure. The uncertainty 
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TABLE II. Light component energy spectrum measured by the ARGO-YBJ experiment by using 
the full 2008-2012 data sample in each energy bin. 


Energy Range 

Energy 

Elux ± total error 

[GeV] 

[GeV] 

[m“^s“^sr“^GeV“^] 

3.55 X 10^ - 5.01 X 10^ 

4.35 X 10^ 

(2.94 ±0.19) X 10“® 

5.01 X 10^ - 7.08 X 10^ 

6.11 X 10^ 

(1.13 ±0.07) X 10-® 

7.08 X 10^ - 1.00 X 10^ 

8.55 X 10^ 

(4.73 ±0.29) X 10-^ 

1.00 X 10"^ - 1.41 X 10^ 

1.21 X 10^ 

(1.94 ±0.12) X 10"^ 

1.41 X 10"^ - 1.99 X 10^ 

1.70 X 10^ 

(7.95 ±0.48) X 10-® 

1.99 X 10^ - 2.82 X 10^ 

2.39 X 10^ 

(3.19 ±0.19) X 10"® 

2.82 X 10^ - 3.98 x 10^ 

3.38 X 10^ 

(1.28 ±0.08) X 10"® 

3.98 X 10^ - 5.62 x lO'^ 

4.77 X 10^ 

(5.07 ±0.31) X 10"^ 

5.62 X 10^ - 7.94 x 10"^ 

6.73 X 10^ 

(2.05 ±0.12) X 10-9 

7.94 X 10^ - 1.12 X 10^ 

9.48 X 10^ 

(8.29 ±0.50) X 10-1° 

1.12 X 10^ - 1.58 X 10^ 

1.33 X 10^ 

(3.40 ±0.21) X 10"^° 

1.58 X 10^ - 2.23 X 10^ 

1.85 X 10® 

(1.43 ±0.11) X 10-^9 

2.23 X 10^ - 3.16 X 10^ 

2.56 X 10® 

(6.24 ±0.49) X 10-11 


on the measured spectrum has been estimated by applying large variations (about 50 %) 
to the fiducial cuts and turns out to be of about 3%. The bins located at the edges of the 
measured energy range are affected by an uncertainty of about ±5%. A variation of the 
quality cuts does not give a significative contribution to the total systematic uncertainty. 


2. Reliability of the detector simulation 

A systematic effect could arise from inaccuracies in the simulation of the detector re¬ 
sponse. The quality of the simulated events has been estimated by comparing the distri¬ 
bution of the observables obtained by applying the same selection criteria to Monte Carlo 
simulations and the data sample collected in each different year. As an example in figure 
the multiplicity distribution obtained from the Monte Carlo events is reported with the mul¬ 
tiplicity distribution of the data. The ratio between the two distributions is also reported 
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FIG. 6. The proton plus helium spectrum measured by the ARGO-YBJ experiment using the full 
2008-2012 data sample. The error bars represent the total uncertainty. Previous measurement 
performed by ARGO-YBJ in a narrower energy range by analyzing a smaller data sample is also 
reported (blue squares) [3j. The green inverted triangles represent the sum of the proton and 
helium spectra measured by the GREAM experiment [2^. The proton (stars) and helium (empty 
stars) spectra measured by the PAMELA experiment |25) are also shown. The light component 
spectra according to the Gaisser-Stanev-Tilav (dashed-dotted line) [2^ and Horandel (dashed line) 
|23j models are also shown. 


showing a good agreement between the two distributions. The contribution to the total 
systematic uncertainty due to the reliability of the detector simulation has been evaluated 
by using the unfolding probabilities and turns out to be about ±6%. 


3. Hadronic interaction models 

In order to estimate effects due to the particular choice of the high energy hadronic 
interaction model in Monte Carlo simulations, a dataset has been generated by using the 
SIBYLL 2.1 model [23 EH]- These data have been compared with the QGSJET dataset 
used in this analysis. In hgurej^the ratio between the multiplicity distributions obtained by 
using QGSJET model and the one obtained by using SIBYLL is reported as a function of 
primary energy. The plot shows that the variation of the multiplicity distributions obtained 
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FIG. 7. Multiplicity distributions of the events selected by the criteria described in section III C 


Values for data and Monte Carlo (black solid line) are reported. The ratio between data and Monte 


Carlo is shown in the lower panel. 





5 5.5 
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FIG. 8. Ratio between the mnltiplicity distributions obtained from QGSJET and SIBYLL based 
Monte Carlo simulations. 


with the two hadronic models is of few percents, giving a negligible effect on the measured 
flux. 


15 









4- Contamination of heavier elements 


A possible systematic effect relies in the contamination of elements heavier than Helium. 
The selection criterion based on the particle density rejects a large fraction of showers pro¬ 
duced by heavy primaries, as shown in hgure The fraction of heavier elements, estimated 
by using the QGSJET-based simulations according to the Horandel model [23], is reduced 
and can be considered as negligible at energies up to 100 TeV. In the lower energy bins the 
contamination is about 1%, whereas in the bins below 100 TeV the contamination does not 
exceed few % and in the higher energy bins it is about 10%. The unfolding procedure has 
been set up in order to take into account the amount of heavier nuclei passing the selection 
criteria. The contribution of this effect is therefore not included in the total systematic 
uncertainty. 

5. Summary of systematic errors 

The total systematic uncertainty was determined by quadratically adding the individual 
contributions. The results are affected by a systematic uncertainty of the order of ±5% in 
the central bins, while the edge bins are affected by a larger systematic uncertainty less than 
± 10 %. 

V. CONCLUSIONS 

The ARGO-YBJ experiment was in operation in its full and stable conhguration for more 
than Eve years: a huge amount of data has been recorded and reconstructed. The peculiar 
characteristics of the detector, like the full-coverage technique, high altitude operation and 
high segmentation and spacetime resolution, allow the detection of showers produced by 
primaries in a wide energy range from a few TeV up to a few hundreds of TeV. Showers 
detected by ARGO-YBJ in the multiplicity range 150 — 50000 strips are mainly produced 
by primaries in the (3 — 300 TeV) energy range. The relation between the shower size 
spectrum and the cosmic ray energy spectrum has been established by using an unfolding 
method based on the Bayes theorem. The unfolding procedure has been performed on the 
data collected during each year and on the full data sample. The resulting energy spectrum 
spans the energy range 3 — 300 TeV, giving a spectral index 7 = —2.64±0.01, which is in very 
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good agreement with the spectral indices obtained by analyzing the sample collected during 
each year, therefore demonstrating the excellent stability of the detector over a long period. 
The resulting spectral indices are also in good agreement with the one obtained by analyzing 
the first data taken with the detector in its full configuration |1]. Special care was devoted 
to the determination of the uncertainties affecting the measured spectrum. The uncertainty 
on the results is due to systematic effects of the order of ±5% in the central energy bins. 
This measurement demonstrates the possibility to explore the cosmic ray properties down to 
the TeV region with a ground-based experiment, giving at present one of the most accurate 
measurement of the cosmic ray proton plus helium energy spectrum in the multi-TeV region. 
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